Andreev reflections in the pseudogap state of cuprate superconductors 
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We propose that, if tlie pseudogap state in the cuprate superconductors can be described in terms 
of the phase-incoherent preformed pairs, there should exist Andreev reflection from these pairs even 
above the superconducting transition temperature, Tc. After giving qualitative arguments for this 
effect, we present more quantitative calculations based on the Bogoliubov-de Gennes equation. 
Experimental observations of the effects of Andreev reflections above Tc — such as an enhanced 
tunneling conductance below the gap along the copper oxide plane — could provide unambiguous 
evidence for the preformed pairs in the pseudogap state. 



Underdoped to slightly overdoped cuprate supercon- 
ductors exhibit a progressive transfer, as the tempera- 
ture is lowered, of the low-energy spectral weight to the 
higher energy region in both the spin and charge chan- 
nels [0j|]. The resulting "pseudogap" begins to appear 
at a cross-over temperature T*, well above the super- 
conducting transition temperature, T^ and the nature of 
this pseudogap state has been the focus of considerable 
recent interest. Many proposals have been presented to 
explain this phenomenon, including various forms of elec- 
tron pairing &fl, antiferromagnetic fluctuations [0,0], 
and spin-charge separation scenarios pj9] , but no consen- 
sus on the origin of the pseudogap has emerged. Among 
the electron pairing scenarios, the phase- incoherent pre- 
formed pair (PIPP), or, normal state pair, scenario sug- 
gested by Emery and Kivelson [0 seems particularly in- 
triguing, since the amplitude of the pairing order param- 
eter is non- vanishing even above Tc. In this scenario, 
the PIPP begin to appear at T* , but do not produce 
superconductivity until Tc because the thermally excited 
phase fluctuations of the pairing order parameter destroy 
phase coherence and hence superconductivity for T > Tc. 
Several earlier experiments, as well as recent teraherz 
spectroscopy results ||l^, seem to support the idea of 
the PIPP, but their existence has not been established 
10,0. In the present paper, we show that a crucial test 
of the PIPP scenario will be the possible Andreev reflec- 
tion (AR) [|l2|,|3| from the pseudogap state: We propose 
that if the pseudogap state contains PIPP, AR will oc- 
cur and will produce, among other effects, an enhanced 
conductance even above Tc that should be readily mea- 
surable in tunneling conductance experiments. 

Two qualitative intuitive arguments suggest that AR 
should occur from a state containing PIPP. First, AR 
(or "retroreflection of a hole") in a normal metal- 
superconductor (NS) junction occurs because an elec- 
tron, which is a well-defined elementary excitation (quasi- 



particle) in the normal metal, is not an elementary ex- 
citation of the superconductor. Therefore, an electron 
from the normal metal, upon entering the superconduc- 
tor, has to "reconstruct" itself as a linear combination 
of the quasiparticles of the superconductor, namely, the 
Bogoliubovons. This requires an electron of momentum 
k and spin a to find a time-reversed mate electron of 
momentum — k and spin —a to form a pair. The hole 
left behind by the mate electron then retraces the path 
of the incoming electron, which is the AR. Naively, the 
same underlying physics should apply to the (assumed) 
preformed pairs in the pseudogap state, because an elec- 
tron is not an elementary excitation of that state either, 
once the pairs are formed, independent of their phase 
coherence. The second heuristic argument involves the 
phase stiffness. The superfluid density in the static, long 
wavelength limit, ps(T), which also determines the su- 
perconductor's phase stiffness, Q{T) — Ti ps{T)/kB, van- 
ishes for T > Tc- In the dynamic Kosterlitz-Thouless- 
Berezinsky (KTB) theory, Ps{i^,T) ^ u)Tc/{1 -f ojTc) > 
as Lu is increased from 0, where Tc is the coherence time 
WM. Precisely this finite frequency phase stiffness, mea- 
sured by a teraherz spectroscopy, was recently reported 
to scale as Qktb — {8/7r)TKTB, in accord with the KTB 
theory of the classical phase fluctuations |10[ . Tktb — Tc 
for two-dimensional (2D) superfluids. Importantly, the 
superfluid density at a finite wavevector (and w = 0), 
Ps{q, T) does not vanish identically above Tc but is given 
by ps{q,T) ^ (<7^+)^, where £,+ is the phase correlation 
length for T > Tktb [jljl ■ AR probes the superfluid den- 
sity at the finite wavenumber q ^ l/^ because it occurs 
within the pairing correlation length, ^, around the in- 
terface between, say, a NS junction. Hence, AR, being a 
proximity effect, is not particularly sensitive to the long 
wavelength physics, which determines Tc, and is expected 
to be present above Tc because p{q, T) ~ (^_|_/^)^ > for 
T* >T > Tc. 



To provide more quantitative calculations of these 
qualitative arguments, we will apply the Bogoliubov-de 
Gennes (BdG) equation [^ with appropriate boundary 
conditions ||l^ to a pseudogap state assumed to be de- 
scribed in terms of the PIPP M. As expected, we find 
that the AR is still present well above Tc and its effects 
may be observed in tunneling experiments. The tunnel- 
ing conductance, G = dl/dV, as a function of the bias 
voltage, V, measured along the in-plane {100} direction 
of d-wave superconductors for a small tunneling barrier 
exhibits a conductance enhancement (CE) inside the su- 
perconducting gap [ p^Jl9| . It is enhanced by the factor 
of 2 in the low T and small tunneling barrier limit with 
a edge at the pairing amplitude [gO|. We will show that, 
if the pseudogap state can be described in terms of the 
PIPP, this CE will also be present above, as well as be- 
low, Tc. 

In a paired state, the quasiparticles (Bogoliubovons) 
can be written in terms of a two-component spinor as 

T^Cx) =1 ) N 1 5 where /(x) and g(x) are, respectively, 

the electron and hole components and obey the BdG 
equation |l^. Inspection of the BdG equation shows 
that ip{x) oscillates on a length scale oi kp^, where kp is 
the Fermi wavenumber, so we introduce the transforma- 

•^W ^ = e^''^-^ ( ''!^M, where w(x) and w(x) 
g(x) J \ v{x) y "- ' ^ ^ 

now vary on a much longer length scale than k^p. For 
kp£^ ^ 1 the BdG equation then becomes ||l^ , po| ^ pi| 



tion 



Eu{-x.) = k_F • Vw(x) -f A(k_F,x)w(x), 

m 

Ev{x) = —kp ■\/v{x) + A*{kp,x)u{x). 



(1) 



We write the pairing order parameter as A(ki?,x) = 
Afc exp[z0(x)]9(x), where A^ = Ad[cos{kxa) — cos(fcj^a)], 
a the lattice constant of a copper oxide plane, and Q{x) 
is the step function, x = {x,y,z) is a 3D vector, where 
we take x to be normal to the interface and z parallel 
to the interface and normal to the copper oxide plane. 
The X > and a; < sides correspond, respectively, 
the paired and normal states. In the PIPP scenario, the 
Cooper pairs are formed above T^ but the phase fluc- 
tuations of the pairing order parameter destroy phase 
coherence and hence superconductivity above Tc. Such 
phase fluctuations are taken to be thermal and static as 
in the previous works y,|2^ . Other kinds of phase fluc- 
tuations such as quantum or temporal ones were argued 
to be less important ^,^ . The thermal fluctuations de- 
stroying superconductivity above Tc are due to a plasma 
of unbound vortices and antivortices, and are described 
by the 2D XY model Q . Each vortex is surrounded by a 
circulating supercurrent which is related with the order 
parameter phase 9{x) by Vs(x) = hWd{x)/{2m), where 
Vs(x) is the local supcrfluid velocity. 
Eq. (|l|) can be solved by taking 



m(x) 
v{x) 



= e 



ik-x 



ve 



.10/2 
-ie/2 



(2) 



to yield a Doppler-shifted local quasiparticle excitation 
spectrum of 



E = kp- V3(x) + Jikp ■ k/m)2 + A2 



(3) 



In writing Eq. (0), we have assumed that the phase 0(x) 
varies slowly on the scale of a single pair size ^ in the sense 
that it is meaningful to specify the energy vs. wavevector 
relation of a quasiparticle at the position x. This condi- 
tion can be written as kp ■ Vs(x) < A^, which is satisfied 
in the pseudogap state as estimated by Franz and Millis 
p2[ , as we will discuss later. The change in the local 
quasiparticle excitation spectrum of Eq. (|3|) will affect 
the spectral propoerties of the superfluid in that a physi- 
cal observable must be averaged over the positions of the 
fluctuating vortices. A similar situation arises, as noted 
by Volovik |23] in the mixed state of a d-wave supercond- 
cutor where the superflow around the field-induced vor- 
tices lead to the residual DOS proposional to ^/H. This 
contributes ~ T\/H to the electronic specific heat, and 
leads to general scaling relations pj] . This suggests that 
the semiclassical approximation can be justified in the 
mixed state. In the present case, we consider a plasma of 
thermally induced vortices instead of a regular Abrikosov 
lattice of field-induced vortices. The essential physics, 
however, remains unaltered because we are interested in 
the q ^ 1/^ scale physics. 

To determine how the AR will affect physical observ- 
ables, we begin with the current, /, which is given by 



I = Im [ei)*V'ip] = Im [e {f*Vf + g*Vg)] 



(4) 



To obtain the quasiparticle wavefunctions -0 of a NS junc- 
tion, we follow Blonder et al. ]l7| and write 



V'A'(x) = 



A 



Jkh -x 



B 



o»ke 



(5) 



for the normal metal, where kg = {—kex, key, kez), and 
je/2 \ / „ pie/2 



V'5(X) = C ( ^"t-.«/2 



e'^-'' + D{ J_-;_,,/2 )e^'"'-" (6) 



for the paired state. Here, kg — (kp + k)kp, k — mE/kp, 
kfi — kp + k', qe = kp + q, and q/i = kp -I- q'. When the 
surface is along the {100} direction, one can show from 
the BdG equation that w+ = u^ — uq and v-^- = v- = vq, 
where 



"0 



vo 




(7) 



Here, rj = kp ■ Vg and A^ = AdCos(2(/)) because we 
consider the {100} surface, where cj) = ta,n~^{kFy/kFx)- 
The energy dispersion relation of Eq. (^) determines only 
the components of the momenta along the lip direction 
of the incoming electrons such that k' • kp/m = —E, 
q • kp/m = \/{E — rj^ — A|,, and q' ■ kp = — q • k^^. 
The components normal to kp are determined by the 
boundary conditions, as can easily be seen. 
Applying the boundary conditions at .x = 0, 

ijNix = 0)-yJs{x = 0) = 0, 
^^^-{x = 0)-^^{x^0)^ 2mVbi>Nix = 0), (8) 



dx 



dx 



where Vj, is the barrier potential energy at the boundary, 
we obtain for {100} tunneling 

cos^ 9 
B^-^Ul- r^) {Z{Z + i cos 61) + {7j/2)^ cos2 0} 



r(^/2)cOs2 



i^o = (1 - r2) [z^ + {jj/2f cos^ e]+{i + ?i) cos^ e, (9) 

where F = vq/uq, with uq and vq given by Eq. (^, 9 
the tunnelig angle, and Z is the dimensionless barrier 
strength given by Z = mVb/kp. In deriving Eq. (|^), 
we have made the approximation rj^ = mvsx/kpx ~ 
r]/{2ep) = rj. The average over the tunneling angle 
was performed with an equal probability appropriate for 
small Z @. Inserting Eq. (|) into Eq. (|), we find that 
the tunneling current is given by I{E) ex 1 + \A\'^ — \B\'^ . 
This current must be averaged over the phase fluctua- 
tions of the order parameter in the pseudogap state as 
I{V) = J dr]P{T])I{V,r]), where P^r]) is the probability 
distribution of 77 given by P{ri) — {S{ri — kp • Vs(x)). The 
angular brackets indicate thermodynamic average over 
the phase fluctuations governed by the 2D XY model 
p2| . In the cumulant expansion, P yields a Gaussian 
distribution of the form 



P{V) 



1 



,-q'l(2-W) 



^/2ttW 



where W is given by 

W w 3.48(aL + aT)A^(r/re). 



(10) 



(11) 



It is crucial to note that ^(77) was evaluated above Tktb 
so that we are in the non-superconducting pseudogap 
state. The a^ comes from the longitudinal fluctuations 
and will be strongly suppressed in a realistic model by 
the Coulomb interaction, while ar comes from the trans- 
verse fluctuations due to vortices p3,pq| . Franz and Mil- 
lis found ut « 0.1 [ p2[ b y fitting Eq. ( [ll| ) to the scan- 
ning tunneling spectra ||26{ . We use this value below. The 
width of the phase fluctuation given by Eq. ( p^ ) and ( |ll| ) 



is VW « 0.6Ad which means that rj < Ad is satisfied in 
the pseudogap state of the cuprates. We, therfore, ex- 
pect that there will exist AR and residual CE, as shown 
explicitly below. If, on the other hand, the phase fluctu- 
ations are too strong in the sense that rj ^ A^, then the 
semiclassical approximation employed to write Eqs. (J^) 
and (P) is not valid, and AR will be washed away. 

From the current we can calculate the tunneling con- 
ductance, G — dl/dV, which after averaging over the 
phase fluctuations, is given by 



1 /•^/2 /-oo 1 

G{V) = - / de dr] , e 



7r/2 



V2^m' 



'''/(2i^)G(F,r,,0), 



GiV,r^,9) 



e r .^df{E-V) 



dE- 



dV 



[1 + \A\^-\B\% (12) 



where A{E,i],9) and B{E,r],9) are, respectively, the co- 
efficients for the Andreev and normal reflections, and 
are given by Eqs. (0) and (g). In Fig. 1, we show 
the normalized tunneling conductance, G{V)/Gn{V) — 
7r(l - Z/Vl + Z^)-^G{V), where G,, is the tunneling 
conductance at the normal state for T > T* averaged 
over the angle 6*, calculated with Eqs. (^) and (p^, as a 
function of the bias voltage. The curves are (from above) 
for T = 1, 1.5, 2, and 3 T^. We take W = 0.3 AI{T/tJ) 
and Ad = 0.04 eV for the underdoped T^ = 84 K BSCCO 
analyzed by Franz and Millis 123], with ep/Ad — 10. Z 
is to be regarded as a fitting parameter to experimental 
data; we have taken Z = 0.1 as a representative value for 
low barrier tunneling experiments like the point contact 
spectroscopies [Q^ ■ The calculated conductance is not 
very sensitive to the precise parameter values. Note that, 
as anticipated, the contribution to the tunneling conduc- 
tance from the AR does not vanish upon averaging over 
the phase fluctuations because we average |Ap and \B\'^. 
In Fig. 1, there is a strong CE due to the non- vanishing 
pairing amplitude, which is clearly seen up to the high- 
est temperature and arises from the AR off the PIPP. At 
T = 3Te, G{0)/Gn{0) ~ 1.46, for the parameters given 
above. Without the normal state pairing, G(0)/G„(0) 
would be equal to 1 above Tc- 

What is the experimental situation, and what further 
improvements in or extensions to the theory should be 
considered? First, experimentally, no CE has yet been 
reported in tunneling experiments above Tc. Our re- 
sults strongly suggest a more systematic search for CE 
for T > Tc in heavily underdoped and high quality clean 
samples. In the heavily underdoped materials, Tc is sub- 
stantially reduced, which implies that the pseudogap phe- 
nomenon can be probed at the relatively low temperature 
where the thermal fluctuations are suppressed. Similarly, 
high quality clean samples will also help to observe the 
CE because the increased electron mean free path will 
enhance the AR |g^. Second, here as in most previ- 
ous calculations, the pairing gap and the wavefunctions 
of corresponding quasiparticles were not determined self- 
consistently, because the previous non-self-consistent so- 
lutions of the BdG equation produced satisfactory results 



for the tunneling conductances |g^ . A main consequence 
of self-consistent solutions is that the pairing amplitude 
varies smoothly from to a bulk value over the corre- 
lation length ^ unlike non-sclf-consistent case where the 
pairing amplitude changes abruptly from to the bulk 
value at the interface. The phase of the pairing order 
parameter, however, remains to be describable in terms 
of the 2D XY model, and will be insensitive to a partic- 
ular way in which the pairing amplitude is varied, and 
hence to the self-consistency. It will neverthless be of 
interest to study the modifications produced by includ- 
ing self-consistency. Finally, our result that a CE can be 
caused by the AR from the PIPP, which are a particle- 
particle (pp) condensate, suggests that one consider pos- 
sible AR from other forms of condensates. AR is novel in 
that it can distinguish pp from particle-hole (ph) conden- 
sates, unlike the most spectroscopic and transport mea- 
surements. For a charge density wave (CDW) or spin 
density wave (SDW) (both ph condensates), we again 
anticipate Andreev-like reflections because an electron 
is not an elementary excitation once the condensate is 
formed. We expect the effects of AR from any ph con- 
densate to be a dip, rather than a peak, around zero 
bias voltage in the tunneling conductance because the 
Andreev-reflected particles are electrons for ph conden- 
sate. We are currently exploring this effect and other 
potential experimental signatures that may distinguish 
among models of the pseudogap state. 

To summarize, we have proposed the Andreev reflec- 
tion as an unambiguous test of the phase incoherent pre- 
formed pair scenario of Emery and Kivelson, which was 
substantiated by the quantitative analyses based on the 
Bogoliubov-de Gennes equation. Experimental observa- 
tion of the Andreev reflection above Tc, such as the en- 
hanced tunneling conductance around zero bias voltage 
along {100} direction, could provide a convincing evi- 
dence for the preformed pairs in the pseudogap state. In 
order to provide more robust argument for our proposal, 
we are currently solving, without assuming the semiclas- 
sical approximation, the Bogoliubov-de Gennes equation 
with the phase dynamics included fully self-consistently. 
This will be reported in a separate paper. 
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Figure Caption 

The normalized tunneling conductance along {100} direction as a function of the bias voltage V in units of the 
pairing amplitude for Z = 0.1. The curves are, from above on the V = axis, for T — I, 1.5, 2, and 3 Tc. The 
conductance enhancement due to the Andreev reflections appears as G{V = 0)/Gn{V = 0) > 1. 



